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Abstract 



We present a new method for the nonhnear approximation of the solution manifolds of 

parameterized nonlinear evolution problems, in particular in hyperbolic regimes with moving 

discontinuities. Given the action of a Lie group on the solution space, the original problem is 

J^ reformulated as a partial differential algebraic equation system by decomposing the solution 

into a group component and a spatial shape component, imposing appropriate algebraic 

constraints on the decomposition. The system is then projected onto a reduced basis space. 

vQ We show that efficient online evaluation of the scheme is possible and study a numerical 

^_^ example showing its strongly iinproved performance in coinparison to a scheme without 

freezing. 
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i/^ 1 Introduction 

.• Reduced basis (RB) methods are effective tools for approximating the solution manifolds of 

f~^ parameterized evolution problems by low-dimensional linear spaces, enabling fast online evaluation 

^f) of the solution for arbitrary parameter values. For linear problems, the POD-Greedy algorithm [5] 

1— I is by now an established reduction approach which has recently been proved to be optimal in the 

^ sense that exponential or algebraic convergence rates of the Kolmogorov n-width are maintained 

•^ by the algorithm [3]. The approach has been further extended to nonlinear settings in [2] based 

rS on empirical interpolation of the nonlinear operators. 

i3 For convection dominated problems with low regularity however, the performance of RB- 
methods is limited by the fact that translation of functions is nonlinear in the translation vector: 
the solution has to be approximated at every location in space it is being transported to, resulting 
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in merely linear decay of the n- width. Our aim is therefore to extend the above approaches by 
additionally allowing transformations of the reduced spaces given by an appropriate group action 
on the solution space. 

Originally developed for the study of relative equilibria of evolution equations pQ [S] , the 
method of freezing allows us to obtain such a decomposition of the solution into a group and 
shape component for arbitrary Lie group actions, provided they satisfy the equivariance condition 
Q. Combining this method with RB-techniques, we obtain a new nonlinear reduction method 
for parameterized nonlinear evolution equations. 

2 The method of freezing 

Assume we are given a parameter dependent nonlinear Cauchy problem of the form 

dtu,,it) + C^{u^{t)) = 0, u,,iO) = UQ (1) 

where M^(t) € F is a function of space for each t g [0, T] and £^ is a partial differential operator 
acting on V. 

Given a Lie group G acting smoothly on V by linear operators, we want to decompose the 
solution Ufj_{t) into a group component ^^(t) G G and a shape component Wp(t) G V such that 

u^{t)=g^{t).v^{t). (2) 

Here the action of G on V is denoted by a lower dot. 

As a guiding example consider the action of the Lie groups G — K*^ on functions v : Mf^ — > M. 
via translation, i.e. {g .v){x) = v{x — g) for g G W'-. In this case, if u^{t) is a solution of M 
drifting along a certain trajectory in space, we want to find a decomposition (pi) such that the 
drift is captured by the evolution oi g^{t) whereas v^(t) becomes as stationary as possible, only 
describing the change of the shape of M^(i) over time. 

Note that, in general, the group action does not have to be induced by a mapping of 
the underlying spatial domain (e.g. shifts or rotations) but can also involve more general 
transformations of V . 

Inserting (pi) into (pi), we formally get 

dtg^{t) . v^{t) + g^{t) . dtv^it) + Cf,{g^{t) . v^itj) = 

which can be rewritten as 

dtv^it) + g^{t)-' . C^{g^{t) . v^(t)) + Q^(t) . v^{t) = 0, Q^{t) = g^{t)-^ ■ dtg^it) (3) 

where 2fi{t) is an element of the Lie algebra LG of G, i.e. the tangential space of G at its neutral 
element Iq- 

Since the decomposition (pi) introduces dim G additional degrees of freedom, the resulting 
system is now underdetermined. It is the main idea of the method of freezing to compensate 
for these degrees of freedom by adding appropriate algebraic constraints <& which force v^ to 
have minimal change over time. These constraints are called phase conditions. Thus, the group 
component of the solution is automatically determined by the phase condition, and no a-priori 
knowledge of the evolution of the solution is necessary. 

If one further assumes that the operator £^t is equivariant under the group action, i.e. 

g-\C^,{g.v)^C^(v) ^g^G^v^V, (4) 



the system decuples into the partial differential algebraic equation (PDAE) of index one 

dtv^it) + Cf,^^^^^){v^{t)) = 0, <P{v^{t),g^it)) = (5) 

with £^[^(u) = Cfi{u) + t) .u and the ordinary differential equation dtgfj.{t) — g^i,{t) ■ Q^{t), called 
the reconstruction equation. The initial conditions are Vfj_(0) = uq and g^(0) = 1g- 

Different choices of phase conditions are possible [1] . We will restrict ourselves here to the 
so-called orthogonality phase condition: assuming that V is equipped with an inner product, 
we require that the evolution of w^ is at each point in time orthogonal to the action of LG, i.e. 
{dtVfj_{t), () . Vfj_{t)) = for all t) € LG. Inserting ^, we obtain 

{C,,iv^it)), f) . v^it)) + (0,,(t) . v^{t), i) . v^{t)) - Vf) e LG. 

After choosing a basis for LG, this leads to a linear dimG x dimG equation system for Q^{t). 
Denoting the basis vectors by Ci, . . . , edimd this system can be written as 



*K(0: 0M W) = (er ■ V^,(t), ts . V^{t)) 



sit) 



{Cf,iv^,{t)), Zr-Vf,it)) 







where g^.s denotes the sth component of fl^ with respect to the chosen basis. 

Depending on the specific problem, many different choices of space-time discretizations of 
the frozen system ([5| are possible. To keep the notation simple, we will restrict ourselves to an 
explicit Euler time-stepping: 

Definition 2.1 Let an H -dimensional space Vh and discrete parameter dependent operators L^, 
^u ti '^'^'^ ^? "" ^H be given, approximating the operators C^, /3p,r) and e^ . ( • )• Since LG acts 
on V by linear operators, we can assume that the §p are linear. Furthermore let P^ : V — > Vh 
be a projection operator. For i^ G N, At = T / K and < fc < K , we define the discrete solutions 
v^ G Vh , 0u G LG of the frozen system (|5L) by the equations wJJ = Ph{uq) and 

with k = 0, . . . ,1 



vt"' = < 



1. Moreover, let 
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A discrete approximation of the solution u^ (t) is then obtained through u*; = 5u ■ "^u i given an 
appropriate discrete approximation of the action of G. 



3 The FrozenRB-scheme 

After freezing and discretizing the original problem (fl]), we will now reduce the resulting system 
^ by projecting it onto an iV-dimensional reduced basis space Vjv C Vf^. To achieve fast online 
evaluation of the reduced scheme, we approximate the nonlinear operators L^, L^p, using the 
method of empirical operator interpolation ^. If (^i, • • 
this method produces indices qi, . . . ,qM and vectors ^i , 



'/i' '^fj..,t) 



, (fiH denotes a basis of the dual of Vh: 
• ■ , Cm, Cf : ■ • • , Cm e Vh such that 
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Given an appropriate projection operator P^ : Vr — > V/v, the fully reduced scheme then reads 



,,fe+i 



= <,^ - AtP^ E™=1 ^,™(L::„. {v%J) ■ ^n 
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(7) 



with v%^^ e Va,, q%^^ e LG, and <_^ == Pn{v°)- 

An appropriate reduced space Vn and the interpolation data Qm, ^m, and ^^ can be constructed 
from solutions of ([6| using reduced basis techniques developed in [5]. Choosing the same 
interpolation points for L^ and Lj^ (^ ensures that expensive evaluations of nonlinear flux functions 
have to be carried out only once for both operators. 

3.1 Offline/online decomposition 

The discrete operators L^, Lj^f, that arise from standard discretizations share the property of 
being local in the sense that for an appropriate basis (pi of Vh each evaluation ipi(hf^{v)) only 
depends on at most C degrees of freedom of v, independently of i, v and H (i/-independent DOF 
dependence [2]). In particular, each of these evaluations can be performed at a speed independent 
of iJ. 

Under this assumption we achieve fast online evaluation of ([7]) using precomputed data as 
follows: 

Proposition 3.1 LeiL^, L^ h have the H -independent DOF dependence property. Then there 
are indices q'l, . . . ,q'j^ with L < 2C'M such that 
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for all V G Vh, m = I, . . . , M . 

//f/'i, . . . , ipN is a basis ofV^, define matrices P, EV, PCL*^'^ and PCR' for 1 < r, s < dimG 
as 

with l<l < L, l<m< M andl<n,n' <N. Let moreover L^, LJ^^ : M^ — > M^^ 



be given as 
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If 'v^ is the coefficient vector of v^ with respect to the reduced basis of Vm , then Vr^ is 
equivalent to 



= v^, -AiP.L^ (EV-v^J 



Jc+l 
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4 Numerical experiment 

As a first test for our new method we consider the Burgers- type problem from 0], i.e. we solve 

dtu^{t) + V ■{hu^[tr)=Q, u^{0)^uo (8) 



for t e [0, 0.3] on the domain O = [0, 2] x [0, 1] identifying {0} x [0, 1] with {2} x [0, 1] and 
[0, 2] x {0} with [0, 2] x {1}. The parameter /x is allowed to vary in the interval [1, 2]. Moreover, 
let b = [ 1 1 ]-^ and Ufi{xi,X2) — 1/2(1 + sin(27rxi) sin(27rX2)). 

Equation ([8| is invariant under the action of the group M^ by translations of fi. Since the 
action of LM^ = M^ is given by negative spatial derivatives, this leads us to the frozen PDAE 






1< r, s < 2 



(9) 



with initial condition v^{0) = uq. The reconstruction equation for (7^(t) is given by dtg^i{t) = Q^{t). 

For the discretization of ^ we choose a finite volume scheme on a 120 x 60 grid and take 100 
steps in time. We use the POD-Greedy and El-Greedy algorithms from [2] on a fixed snapshot 
set for the generation of the reduced basis and the operator interpolation data, choosing for 
simplicity a common interpolation basis ^m = £,m fo^' the operators. 

Figure [l] shows that the FrozenRB-scheme provides already for 20 basis vectors a good 
approximation of v^ (right column) whereas, for the same amount of basis vectors, the solutions of 
the same scheme without freezing are strongly deformed (center column). In Figure^ we compare 
the convergence of the schemes to the corresponding detailed approximations: the FrozenRB- 
scheme improves the approximation error for equal basis sizes by two orders of magnitude. 
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Figure 1: Numerical solutions of the Burgers problem (Is]) for // = 1.5; left column: detailed 
approximation; center column: reduced approximation; right column: reduced approximation of 
the frozen solution w^ (20 basis vectors, 38 interpolation points). 
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Figure 2: Error of the FrozenRB-approximation of (|9]) in comparison with the same reduction 
method without freezing. The number of interpolation points M is given by Af = 1.8 • N. 
The error is estimated as the maximum approximation error over a set of 100 randomly chosen 
parameters. 



